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Abstract 

This paper deals with existence of strange stars in / (T) modified 
gravity. For this purpose, we have taken the diagonal tetrad field of 
static spacetime with charged anisotropic fluid and MIT bag model, 
which provide the linear relation between radial pressure and density 
of the matter. Further, the analysis of the resulting equations have 
been done by assuming the parametric form of the metric functions in 
term of the radial profiles with some unknown constant (introduced 
by Krori and Barua). By the matching of two metrices, unknown 
constant of the metric functions appear in terms of mass, radius and 
charge of the stars, the observed values of these quantities have been 
used for the detail analysis of the the derived model. We have dis¬ 
cuss the regularity, anisotropy, energy conditions, stability and surface 
redshift of the model. 
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1 Introduction 


Recently, the modified theories of gravitation have attracted the attention 
of a huge community of theoretical physicists due to their implication as an 
alternative framework for the successful proposal of dark energy. This modi¬ 
fication is carried out in several ways by modifying the Lagrangian of General 
Relativity (GR), for examples Ricci scalar R in the Einstein-Hilbert action 
must be replaced by / {R) (see Durrer and Maartens 2010; DE Felice and 
Tsujikawa 20210). In GR and modified theories of gravity, the gravitational 
effects could be described in term of contents of matter and curvature of 
spacetime. However, one can formulate an analogous theory of GR, which 
can be demonstrated in term of torsion instead of curvature. Originally, Ein¬ 
stein formulated this theory to unify the basic forces due to Goulomb field 
and gravitational field Einstein (1928). This problem of unification remains 
unsolved in the new theory, but after many years the theory becomes much 
popular as an alternative theory of GR , this theory is known as teleparallel 
equivalent of GR (TEGR) Mdller (1961). The basic concept to formulate this 
theory is to take a such manifold that preserve torsion along with curvature. 
The complete Riemann curvature tensor can be considered as zero, hence 
in this way gravitational field can be described by torsion. The appropriate 
technique is to utilize the tetrad fields Weitzenbdck (1923). 

During the current years, the researchers have carried out the study of 
TEGR with some modification to explore some theoretical as well exper¬ 
imental problems of modern cosmology (Ferraro and Fiorini 2007, 2008). 
The comparison of / {R) gravity with /(T) gravity implies that later one is 
more convenient as in /(T) gravity, equations of motions are of second order 
while for metric / {R) gravity, these equations are fourth order. Recently, 
a lot of work in / (T) gravity is devoted to investigate the observed data of 
the theoretical cosmology or to determine the constraint on the parameter 
of theory, so that the theoretical results may fit with the data obtained by 
some standard source like WMAP etc. Further, many researchers have inves¬ 
tigated that a particular class of / (T) gravity models can explain some basic 
problems of cosmology like inflation and accelerated expansion (Bengochea 
and Ferraro 2009 Wu and Yu 2010). 

In addition to study the cosmology, f{T) theory has been used to exam¬ 
ined the effects of f{T) models on the existence of Black hole and compact 
objects in 3-dimensions. On the basis of this fact, the BTZ black hole model 
has been derived in f{T) theory of gravity. (Zheng and Huang 2011). Also, it 
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has been proved that the first of black hole thermodynamics in f{T) gravity 
becomes invalid, due to the break down of Lorentz invariance (Bengochea 
2011). Further, some static vacuum solutions with charged source have been 
investigated in f(T) theory, such solutions leads to charged black hole so¬ 
lutions in f(T) gravity (Li et al.2011). A large class of static perfect fluid 
solutions exhibit the existence of relativistic stars in /(T) gravity ( Dent et 
ah2011, Rong-Xin et ah2011; Wang et ah 2011, Boehmer et ah 2011, Daouda 
2011 ). 

It has been the subject of great interest to study the models of anisotropic 
stars during the last decades in GR as well in modihed theories of gravity 
(Herrera and Santos 1997, Chan et al.l993, Di Brisco 1997, Herrera and Bar¬ 
reto (2013), Abbas (2014a), (2014b), (2014c), (2014d), Abbas and Sabiullah 
(2014)). Mak and Harko (2004) present a class of regular and analytic so¬ 
lution of held equation with anisotropic source. Also, Chaisi and Maharaj 
(2005) have discussed the nature of anisotropic matter source, after estab¬ 
lishing an algorithm. By utilizing the Chaplygin gas equation of state (EOS), 
Rahaman et ah (2012) generalized the Krori-Barua (KB) (1975) solutions 
to the matter source with electric held. It has become a scientihc tool to 
study the models the compact stars with Krori-Barua metric (Kalam et al. 
2012, Kalam et al. 2013). Hossein et al. (2012) discussed the modelling of 
the compact star model by taking cosmological constant as radial coordinate 
dependent. Bhar et al. (2015) and Maurya et al.(2015) discuss the pos¬ 
sibility for the existence of higher dimensional and charged compact stars, 
respectively. 

In present year, Abbas and his collaborators ( Abbas et al.2014, Abbas et 
al.2015a, 2015b, 2015c, 2015d ) have studied the models of anisotropic com¬ 
pact stars in GR, f{R), f{G) and f(T) theories in diagonal tetrad case using 
Krori and Barua (KB) metric approach. Das et al.(2015) have examined 
the existence of compact stars in /(T) gravity using the conformal motion 
technique and have discussed properties of the resulting models. The main 
objective of the present paper is to study the models of the anisotropic com¬ 
pact stars in the context of f{T) gravity using diagonal tetrad in the presence 
of electric held and MIT Bag Model. We have examined the anisotropic be¬ 
havior, regularity at the center as well as stability of these models. Finally, 
the surface redshift has been calculated. All these properties of the models 
have been discussed by using the observational data of the compact stars. 
The plan of the paper is the following: In the next section, we present the 
review of f{T) gravity and equations of motion with charged anisotropic. 
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Section 3 deals with charged anisotropic source and the equation of motion 
in f{T) gravity with diagonal tetrad. Section 4 investigates with the physical 
analysis of the proposed models. In the last section, we conclude the results 
of the paper. 


2 f(T) Theory of Gravity 

Currently, Teleparallel theory of gravity is an equivalent theory of gravitation 
to the GR (Zhang et al.2011, Li et al. 2011). In this section, we shall 
introduce the basic concepts of the /(T) gravity theory. To this end, we 
dehne notation of the Latin indices for the tetrad held and Greek indices for 
spacetime coordinates. The line element of the manifold is given by 

ds'^ = QfMudx^dx’'. ( 1 ) 

This spacetime can be transformed to a Minkowskian form by the matrix 
transformation as follows: 


dS^ = g^^dx^dx'' = Vij9^e\ ( 2 ) 

9^ = e^^dx^, (3) 

where riij = diag\i, —1, —1, —1] and efej' = 5^ or efe^ = 51- 
The root of the metric determinant is given by = det[e^^] = e. 

We consider the manifold for which the Riemann tensor is zero and non¬ 
zero torsion terms exist, the Weitzenbock’s connection components can be 
dehned as follows: 


T 


a 

fiiy 



(4) 


We dehne the torsion and the contorsion tensor as follows: 


rpCX _ ■pQ; 


- a^e;). 


and the components tensor are dehned as 


sr = §(AT + 




(5) 

( 6 ) 

(7) 
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one can write the torsion scalar as 


T = (8) 

Now, similarly to the f{R) gravity, one dehnes the action of f{T) theory as 
follows: 


5'[e^, ^a] = j d*xe[-^f{T) + CMatter{^ a)], (9) 

where we used G = c = 1 and is matter helds. Now the variation of 
above action yields the following form of equations of motion (Wu and Yu 
2010 and Zheng and Huang 2011) 

SyapT/rr + + T," = 4^ (7J + E;) , (10) 

where 7)]^ and i?)) are energy momentum tensors of ordinary matter and 
electromagnetic field. 

The ordinary matter is an anisotropic fluid for which the energy-momentum 
tensor is given by 

V = (d + Pt)u^,u'' - pth)) + {pr - Pt)Vf,v'', (11) 

where is the four-velocity, radial four vector, p the energy density, 
Pr is the radial pressure and pt is transverse pressure. Further, the energy 
momentum tensor for electromagnetic held is given by 

( 12 ) 

where is the Maxwell held tensor dehned as 

(13) 


and is the four potential. 


3 Equation of Motion in f{T) Gravity 

We assume that the interior spacetime of a strange star is describe by the 
KB (1975) metric 
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+ r‘^{d9‘^ + sin^dcl?) 


(14) 


In order to re-write above line element, into the invariant form under the 
Lorentz transformations, we dehne the tetrad matrix as 

• Q.(r) b{r) 

= diag[e~ ,e~ ,r,rsin{6)]. (15) 

Further, one can obtain e = def[e^] = r‘^sin{9). For the charged fluid 

source with density p(r) radial pressure Pr{r), tangential pressure pt{r). The 
Einstein-Maxwell (EM) equations take the form ( with geometrized units 
G' = c = l). 


T(r) 

2e-\ , 1, 

— ^ + ^)’ 

(16) 

T'(r) 

2e"\ // 1 E, 

— (a -t- T{b + )), 

iY> qn ^ iY> 

(17) 


where the prime ' denotes the derivative with respect to the radial coordinate 
r. One can now re-write the equations of motion for anisotropic fluid as 


dvrp -|- 
Anpr — E'^ 


A'Kpt + E'^ 
cot9 , 

E{r) 






It 

2 


/ 

4’ 




, o 1 , , / 



/t 

2 


= 0 , 



(18) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 


where q (r) is the total charge within a sphere of radius r. According 
to the MIT bag model, we take the simple form of the strange matter EoS 
(Rahman et ah(2012)) 


Pr 


^{p-ABg). 


(23) 
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From Eq. fl2ip . we get the following linear form of f{T) function 

f{T) = f3T + f3,, (24) 

where j3 and /3i are integration constants, for simplicity we assume /3i = 0. 
We parameterize the metric as follows (Krori and Barua 1975) 

h{r) = Ar'^, a{r) = Br"^ + C, (25) 

where A, B and C are are some arbitrary constants which will be determined 
later using some physical conditions. With the choice of the EoS, we have a 
system of hve independent equations with Eve unknown parameters namely, 
p, Pr,Pt, E{r) and (j{r). From Eqs. ffT8|l - fl25|) . we obtain 


P = 

Pr = 
Pt = 

= 


3{A + B) ^^^2 
(71 + S)^_^,2^_ 


-B — -A + B^r^ — ABr^ H—- 


Ibvr 

1 

Stt 

+ Bg, 


The charge density is obtained as 




/5- 




^2 _ 

a = [2 - (.1 - 35) r^] + 


27rr 


Sttvi/j 




ATTr^ijj 


where 


^ = 


^(A - 35 - + ^ - SnBg 


The amount of net charge inside a sphere of radius r is given by 


q = r 




(26) 

(27) 

(28) 
(29) 


(30) 

(31) 


(32) 
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4 Physical Analysis of Model 


4.1 Central Regularity of Model Parameters 


It well know fact that the KB metric preserve no singularity in its metric 
functions even at the center at r = 0. Therefore, it becomes necessary to 
impose the constraints on the constants appearing in the metric functions, so 
that physical parameters of the model remain well behaved at all the inner 
points strange stars. For the regularity at the center (r = 0), we obtain the 
central density in the form 




(33) 


Further, the regularity of electric held requires that it must vanish at the 
center i.e.. 


which implies that 


= 0) = i;(A -B)- SjtB, = 0, 




(34) 


(35) 


^ IGtt 

The pressures and density should be decreasing function of r. In our model, 
radial variation of is obtained as 


dpr 

dr 


(A + B) 

Stt 


< 0 , 


(36) 


At r = 0, ^ = 0 and < 0. This implies that Pr will decrease radially 
outward. The radial variation of matter energy density is obtained as 


(k = (37) 

dr Stt 

which also shows that at r = 0,^ = 0 and ^ + S) < 0. Hence 

density is decreasing function of r, provided that /9 > 0, as we have chosen it 
in our discussion. Hence, it has been shown graphically that pressures and 
energy density decreases as shown in hgures 1-3. The anisotropic stress is 
obtained as 

A = pt-pr = 2Bg-^ 

+ — {3B -2A + B^r^ - ABr^ + 

Stt H 
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(38) 




Figure 1: This figure represents the variation of p versus r{km) for a strange 
star of radius 8.26 km with /3 = 1, /5 = 2 and (5 = 3. The corresponding 
values of A, B and Bg have been used from table 2. All the graphs have 
been plotted for table 2, it will not be mentioned again. 



Figure 2: This figure represents the variation of pr versus r{km) for a strange 
star of radius 8.26 km with /5 = 1, /5 = 2 and (3 = 3. 


The anisotropic force of the gravitating system will be will be outwardly 
directed when pt > Pr i.e., A > 0, while inwardly directed when pt < Pr i.e., 
A < 0. It is clear from the figure 4 that there exist repulsive gravitational 
force for 1 < /? >< 2 as (A > 0), while for (3 >3 force becomes attractive. In 
the case of repulsive force more massive distribution of matter would exist, 
which is the result of disturbance of equilibrium. Figure 5 implies that E'^ is 
increasing function of radial coordinate for all the values of (3. 



\ 

\ 

\ 


Figure 3: This figure represents the variation of pt versus r{km) for a strange 
star of radius 8.26 km with (3 = 1, (3 = 2 and (3 = 3. 
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4.2 Matching Conditions 

The exterior region of the star is taken as Reissner-Nordstrom metric given 
by 

ds^ = -(l-‘^ + 9l'\dt^+(l- — + ^] + r\de^ + sin^Od^'^) 

\ j V j 

(39) 

where Q is the total charge within boundary r = R. Continuity of the metric 
coefficients gtt, grT and across the boundary surface r = R between the 
interior and the exterior regions of the star yields the following results: 


2M _ bb^+c 

R^ R^ 

3 ‘2M 

R^R^ 

^ — no BB?+C 

By solving the above three equations, we get 


^ = -Jpln 


2M g2 


(40) 

(41) 

(42) 


(43) 


1 

'M 



2M g2' 


R2 

R2_ 


1 

1- 

1 

4 

1 

bO 1 
1_ 


C 


In 


2M g2 


M _ 

R _ R? 


1 2M I 


(44) 

(45) 


For the given values of the parameters M, R and Q the values of A and B is 
are given in table 2. 


4.3 Energy conditions 

The anisotropic charged fluid considered as strange matter will satisfy the 
null energy condition (NEC), weak energy condition (WEC) and strong en¬ 
ergy condition (SEC) if the following inequalities hold simultaneously at all 
points for the given radius of the star: 
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Figure 4: This figure represents the variation of A versus r{km) for a strange 
star of radius 8.26 km with /5 = 1, /5 = 2 and (3 = 3. 



Figure 5: This figure represents the variation of versus r{km) for a strange 
star of radius 8.26 km with /5 = 1, /5 = 2 and (3 = 3. 



Figure 6: This figure represents the variation of p + ^ versus r{km) for a 
strange star of radius 8.26 km with f3 = 1, (3 = 2 and (3 = 3. 



Figure 7: This figure represents the variation of p + Pr versus r{km) for a 
strange star of radius 8.26 km with (3 = 1, (3 = 2 and (3 = 3. 
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Case 

M 

R{km) 

M 

R 

W 

1 

IAMq 

8.26 

0.25 

0.004 

2 

IAMq 

6.88 

0.30 

0.027 

3 

IAMq 

5.90 

0.35 

0.061 

4 

IAMq 

5.16 

0.40 

0.105 


Table 1: Values of M, R and Q defined in Rahman et al.(2012) 


Case 

A {km 

B {km 

Bg {km 

1 

0.0102 

0.0073 

0.0001732 

2 

0.01798 

0.01351 

0.000267 

3 

0.0292 

0.0231 

0.0003643 

4 

0.044 

0.037 

0.000418 


Table 2: Values of the model parameters A, B and Bg for given Masses and 
Radii of Stars 


ta2 


P + Pr > 0, 


P + Pt + ^ >0, 
Air 


(46) 

(47) 

(48) 


£2 £2 £2 


p+p,+ — 

4n 


p+p ,+— 

4n 

p+p,+ - 

4n 


0.0014 



0.0028 


0.0042 



0.0013 



0.0027 


0.0041 

" - ^ 


0.0012 


x, 

' \ 

0.0026 


0.0040 


X 

0.0011 


\ 

0.0025 


“V 0.0039 



0.0010 



0.0024 


0.0038 


X 

0.0009 


\ 

0.0023 


0.0037 


X 



. 





X 



4 6 ^ . 

' 

1 4 6 K ' 

2 4 6 gs. 


Figure 8: This hgure represents the variation of P + Pr + ^ versus r{km) for 
a strange star of radius 8.26 km with = 1, /5 = 2 and /3 = 3. 
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Case 

p (r = 0) {gm cm 

p{r = R) {gm cm 

Pr {dyne cm 

1 

1.643x10^^ 

9.362x10^^ 

2.123x10^^ 

2 

2.895x10^^ 

1.443x10^^ 

4.361x10^^ 

3 

4.703 xlO^^ 

2.015x10^^ 

8.204x10^^ 

4 

7.087x10^^ 

2.58x10^'^ 

14.48x10^^ 


Table 3: Values of density and pressure at center and surface of star for /5 = 1 


Case 

p (r = 0) {gm cm 

p{r = R) {gm cm 

Pr- {dyne cm 

1 

1.643 xlO^^ 

9.362x10^4 

2.123x1055 

2 

2.895x10^^ 

1.443x1045 

4.361x1055 

3 

4.703 xlO^^ 

2.015x1045 

8.204x1055 

4 

7.087x10^^ 

2.58x1045 

14.48x1055 


Table 4: Values of density and pressure at center and surface of star for (5 = 2 


p + Pr + 2pt + — >0, (49) 

47r 

Employing these energy conditions at the centre (r = 0), we get the following 
bounds on the constants A and B\ 


£j2 

{i)NEC ■p + —>0^A>0. 


£j2 

{ii)WEC :p + pr>0^A + B>0,p + pt + —>0^A + B>0 



Figure 9: This figure represents the variation of p + Pr + 2pt + ^ versus 
r{km) for a strange star of radius 8.26 km with (3 = 1, (3 = 2 and (3 = 3. 
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Case 

p (r = 0) {gm cm ^) 

p{r = R) {gm cm ^) 

p,, {dyne cm 

1 

1.643x10^^ 

9.362x10^^ 

2.123x10^^ 

2 

2.895x10^^ 

1.443 xlO^^ 

4.361x10-^^ 

3 

4.703 xlO^^ 

2.015x10^^ 

8.204x10^^ 

4 

7.087x10^^ 

2.58x10^^ 

14.48x10-^^ 


Table 5: Values of density and pressure at center and surface of star for /3 = 3 

j^2 

{iii)SEC : p + Pr + ‘^Pt H-> 0 ^ -B > 0. 

4:71 

Since > 0, so condition (i) is satisfied. The weak and strong energy 
conditions (ii) and (iii) will be satisfied if I? > 0. With the set of values given 
in table 1 and 2, we have shown in Figures 6-9 that the energy conditions 
are valid through out the interior of the given star. 

4.4 Stability 

In order to discuss the stability of our constructed models, we follow the ap¬ 
proach based on Herrera (1992) cracking (or overturning) concept. According 
to this technique the squares of the radial and tangential sonnd speeds must 
lies in the interval [0,1]. Farther, cracking concepts predicts that the region 
for which the radial speed of sonnd is greater than that of transverse speed 
, then snch region is potentially stable region. It is clear that, for no crack¬ 
ing (stability), the difference of two sonnd speeds, i.e.,nf^ — should attain 
the same sign everywhere inside the anisotropic matter distribntion. In our 
model, we have 


where. 


dpr 1 

dp 3’ 


2 dpt a + j 


a = e-^^"/3 


(^B^r — 2ABr — 



1 


(50) 

(51) 


(52) 
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Figure 10: This figure represents the variation of v‘1^ — v‘1^ versus r{km) for 
a strange star of radius 8.26 km with = 1, /5 = 2 and /3 = 3. 


7 = -2Are-^^ 13 


^B — -A + — ABr^ H—r 


^2 2 r 

For causality condition to be satisfied, we must have 

a + 7 


0 < 


< 1 


(53) 


(54) 


-3(A + B)rAf3e-^''^ 

For the assumed set of values, figure 10 shows that the above condition 
is satisfied. Hence, our proposed model preserve stability in f(T) gravity . 


4.5 Surface Redshift 

We define the compactification factor as 




R 


The surface red-shift {Zg) corresponding to the above compactness {u) is 
obtained as 

Z, = (l-/3(l-e-^^")Y" -1 


The maximum surface red-shift, in this set up, for a strange star of mass 
1.4M0 and radius 8.26km turns out to be Zg = 1.4 for /5 = 1, while Zg < 1.4 
for /5 > 1 as shown in Figure 11. 


5 Discussion 

During the last few years, theoretical physicists have shown their attraction 
to study gravitational field as the effect of torsion instead of curvature of 
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R 


Figure 11: This figure represents the variation of redshift Zg versus r{km) 
for a strange star of radius 8.26 km with /5 = 1, /? = 2 and /9 = 3. 


the underlying geometry. The torsion theory was originally derived in par¬ 
allel lines to general theory of relativity which describes gravity in term of 
curvature. Hence, this idea was promoted as teleparallel equivalence of GR. 
The black hole (BH) and neutron stars are the highly dense astronomical 
objects which store information about the entropy on the BH horizon. In 
the modified TEGR, f{T) gravity posses many significant feature as com¬ 
pared to GR. Recent observations from solar system orbital motions in order 
to constrain f{T) gravity have been made and interesting results have been 
found. The general constraints for the existence/non-existence of compact 
stars have been studied in /(T) gravity (Boehmer et ah 2011). It is inter¬ 
esting to drive the models of the compact stars in /(T) gravity using the 
diagonal tetrad field with charged anisotropic fluid and MIT Bag model. 

This paper deals with construction of analytic models of compact stars in 
f{T) gravity with charged anisotropic matter and MIT Bag EoS. The inner 
region of the stars have been taken as static charged anisotropic spherical 
source. By using the diagonal tetrad field, the equations of motion have been 
formulated to complete the discussion. One of the field equations implies that 
unknown function f{T) appears as a linear function of T i.e., f{T) = /3T+I3i, 
where jS and /5i are the constant of integration. Using this form of f{T) with 
KB metric functions, we have determined explicitly the matter components 
and electric field intensity. The anisotropy, regularity and energy conditions 
have been discussed in detail. The observed values of masses, radii and charge 
of compact stars have been used to calculate the unknown constants of KB 
metric. The first and second derivatives of density and pressures, implies the 
maximality of these quantities at the center, and these have decreasing radial 
profile. The values of density and pressure at r = 0 and r = R are given in 
tables 3-5. 

We have found that (see Figure 4) that there exist repulsive gravitational 
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force for 1 < /3 < 2 as (A > 0), while for (3 > 3 force becomes attractive. The 
hrst case regarding to repulsive gravitational force leads to the construction 
of more massive object, while second case gives the formation of smaller 
objects. The Herrera’s cracking concept states that the region for which the 
radial speed of sound is greater than that of transverse speed, then such 
region is potentially stable region. It is clear that for stability, the difference 
of two sound speeds, i.e.,Vg^ — should attain the same sign everywhere 
inside the anisotropic matter distribution. In our discussion, the compact 
stars remains stable in /(T) gravity even in the presence of electromagnetic 
field. The maximum surface redshift, for a strange star of mass IAMq and 
radius 8.26km in /(T) gravity turns out to be Zg = 1.4 for (3 = 1, while 
Zg < 1.4 for /3 > 1 as shown in Figure 11. 
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